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Non-thermal processes in the cluster of galaxies Abell 3376 
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ABSTRACT 

We model the high-energy emission that results from the interaction of relativistic 
particles with photons and matter in the cluster of galaxies Abell 3376. The presence 
of relativistic particles is inferred from the recently found radio relics in this clus- 
ter, being one of the most prominent examples of double opposite, giant ringlike radio 
structures. Assuming that diffusive shock acceleration takes place in the cluster regions 
where radio relics are observed, we calculate the spectral energy distribution result- 
ing from the most relevant non-thermal processes, which are synchrotron radiation, 
inverse Compton scattering, relativistic Brcmsstrahlung, and inelastic proton-proton 
collisions. In the context of our model, the major radiative component at high energies 
is inverse Compton scattering, which could reach luminosities L '--^ 9 x 10^^ erg s~^ in 
the energy range between ~ 1 MeV and 10 TeV. Hadronic interactions would yield a 
minor contribution to the overall non-thermal emission, but would dominate at ultra- 
high energies. The cluster Abell 3376 might be detectable at gamma-rays by HESS, 
GLAST satellite and future planned Cherenkov arrays. 
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1 INTRODUCTION 

Clusters of galaxies are characterized by thermal emission in 
the 2 to 10 keV soft X-ray region (Lx ~ lO"^ - lO"^ erg s"^) 
as a result of Bremsstrahlung radiation of the hot intra- 
cluster gas {kT ~ 5 — 8 keV). However, the detection of 
radia tion at radio wavelengths IjFeretti. Burigana fc En/jlinl 
|2004| '1. as well as in the ex treme ultraviolet (0.07 — 0.4 keV; 
EUV) teowyor et al."2004') and at hard X-rays (20-80 keV; 
HXR) (JFusco-Fcmiano ct al. 2004), indicates the presence of 
non-thermal activity in these systems. There is general con- 
sensus that radio emission is produced by synchrotron radi- 
ation originated from the interaction of rel ativistic electrons 
with the cluster magnetic field {B ~ /^G) (jGovoni fc Ferettil 
[2004). On the other hand, EUV radiation might be pro- 
duced either by a cooler therm al component {kT ~ 2 keV) 
iLieu et al.lll996l : lMittaj[l993 ) or by inverse Compton (IC) 
scattering of cosmic microwave background (CMB) photons 
by the same populatio n of r elativistic electrons res ponsible 
for th e radio emission (|En/31in fc Bicrmann 1998 ; Lieu et al.l 
1 19991 ). The latter process might also explain the HXR emis- 
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sion (|Fusco-Femiano et al.ll 19991 ) . although this possibility is 
still matter of debate. 

Regarding the radio emission, whose non-thermal origin 
is firmly established, clusters of galaxies present two different 
sources of diffuse large-scale synchrotron emission, known as 
'radio halos' and 'radio relics'. The former are located at the 
centre of clusters and are characterized by unpolarized radio 
emission with structures that roughly resemble that shown 
by the X-ray emission, whereas the latter are polarized radio 
sources more irregularly shaped found at the periphery of 
clusters (e.g. iFeretti fc Giovanninill2008l ). Radio relics have 
been observed in several clusters. The most extended and 
powerful sources of this class have been de tected in clusters 
wit h central radio halos, such as Coma (IGiovanriini et al.l 
1991) and Ab ell clu sters A216 3 iFeretti et al.ll200ll), A 2255 
(.Fcretti et"ai] ll997D A22 56 l|Rottgering et al.l 1 19941 ) and 
A2744 l|Govoni et alll200ll ). Only few clusters present dou- 
ble opposite relics, be ing the most promine nt examples 
those found i n A36 67 (|Rottgering et all 1 19971 ) . and A3376 
iBagchi et al.l (|2006l ). The large radio flux recently detected 
from the relics of A3376 is an evidence of that relativistic 
leptons are present in those regions. For this reason, in the 
present work we will be concerned with the non-thermal 
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processes taking place in the outskirts of the cluster Abell 
3376 where the two giant and almost symmetric radio relics 
are located. Nevertheless, the model that we will present in 
this paper is also valid for clusters which present only one 
relic. In the case of Abell 3376 the radio flux is more than 
one order the magnitude larger than in other radio emitting 
clusters, like Abell 521 and AbeU 2163. 

There is increasing evidence that the radio emitting 
particles in relics are accelerated by accretion and merger 
shocks generated in the intracluster medium (ICM) dur- 
ing cosmological large-scale structure formation. Because 
of the electron short radiative lifetimes, radio emission can 
be efficiently produced close to the location of the shock 
waves. Thus, the origin of radio relics can be explained if 
they ar e considered to trace the position of t hese very large 
shocks l|En/?lin et al.lll998l : iHoeft et al.ll2004l '). According to 
the standard diffusive shock acceleration theory, shock waves 
in the presence of even modest and turbulent magneti c fields 
are s ites of efhcient acceleration of charged particles (|Drurvl 
119831 ). Another possibility may be adiabatic compression in 
these environmental shock waves (JEn/jlin fc Gopal-Krishnal 
I2OOII ). 

Based on the inferred population of high energy parti- 
cles involved in the previously mentioned non-thermal pro- 
cesses, it is natural to expect gamma-ray emission from 
galaxy clusters. This emission could be generated from 
neutral pion decay in cosmic ray coU isions in the ICM 
ijVolk. Aharonian fc Breitschwerdtlll996n . or by IC scatter- 
ing that involves CM B photons and relativistic electrons 
IjAtovan fc Volkll200Cl ). Hydrodynamical cosmological sim- 
ulations are a useful tool to estimate the properties of the 
cosmic ray populations in galaxy clusters and their effects 
on thermal cluster observables, as well as to predict the radi- 
ation emitted by these large virializcd objects at 7-ray ener- 
gies (Kcshct 2003, Pfrommcr ct al. 2007a, P frommer et al.l 
20081). A different numerica l approach has been used by 
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evolution of particle and photon spectra resulting from non- 
thermal processes at the shock fronts formed in merging 
clusters of galaxies. All these works point to clusters in gen- 
eral as sources of 7-ray emission observable with GLAST 
satellite. 

These theoretical results are consistent with the fact 
that no 7-ray emission has been detected so far with the 
current observational facilities, as indicated by the lack of 
correlations between unresolved EGR ET 7-ray sources an d 
nearby X-ray bright galaxy clusters (IReimer et al.l 120031 ) . 
Only marginal evidence for emission from A1758 within the 
location erro r contours of t he so urce 3 EG J1337-f 5029 was 
reported by iFegan et all (120051 ). Perseus and A bell 2029 
galaxy clusters were observed by IPerkina (|2006l ) with the 
Whipple 10 m Cherenkov telescope. They find no evidence 
of point source or extended 7-ray emission in the TeV energy 
range. In addition, rece nt observations o f the c lusters Coma 
and Abell 496 made bv lDomainko et al.l (|2007|) with the ar- 
ray of Chrerenkov telescopes HESS have not detected signif- 
icant signals of 7-ray emission in exposure times of ~ 10 — 20 
hours. 

In the present work, we estimate the high energy emis- 
sion of the recently detected radio relics in the nearby cluster 
Abell 3376. Their radio power at 1.4 GHz, together with the 
assumption of an equipartition magnetic field, constrain ef- 



fectively our model which is specific of this particular cluster. 
We assume a diffusive shock acceleration mechanism, and 
include different non-thermal processes that may be acting 
in the radio structures. We predict that the 7-ray photons 
created in this cluster should be detected in the near future 
by GLAST, as well as by the new generation of Cherenkov 
arrays in the southern hemisphere, such as HESS II. 

The paper is organized as follows. Section [2] describes 
the main features of the cluster A3376. Section |3] presents 
the acceleration and loss mechanisms that affect the con- 
tent of relativistic particles in the ICM. The estimates of 
production of 7-rays and lower energy radiation are given 
in Section m Finally, we summarize our results and present 
our concluding remarks in Section [5] 



2 THE CLUSTER ABELL 3376 

The rich cluster of galaxies Abell 3376 has been detected 
by ROSAT and J^CMM-Newton through its X-ray emis- 
sion revealing strong evidence for merger activity of sub- 
clusters. It has a X-ray luminosity Lx(0.1 — 2.4 keV) ~ 
1.22 X 10*** h"'^ ergs^^. This power is produced by ther- 
mal Bremsstrahlung radiation from the hot ICM, which 
has an overall temperature Tx ~ 5.8 x lO'' K (~ 5 keV). 
This cluster is located in the southern hemisphere (a = 
06''0™43", S = -40°03') at a redshift z ^ 0.046, which 
corresponds to a distance of d ~ 197 Mpc for a standard 
A cold dark matter (ACDM) cosmology (r2ni=0.3, nA=0.7, 
Ho = 100/ikms"^ Mpc"\ with h = 0.7). The cluster mass 
has been estimated by applying the virial theorem to the 
cluster member gal axies, assuming th at mass follows the 
galaxy distribution IjGirardi et. allll998r ). giving a virial mass 
Mvir ~ 3.64 X 10^"* h~^ Mq. The corresponding virial radius 
is i?vir ~ 0.98 /i"^ Mpc. 

Radio observations made bv lBagchi et al.l l|2006l ) with 
the Very Large Array (VLA) instrument show the pres- 
ence of two giant, ring-shaped structures in the periph- 
ery of the cluster at a distance of ~ 0.7/i~^ Mpc from 
its centre. The radio fiux detected at i^ = 1.4 GHz is 
Fv = 302 mjy, which corresponds to a radio luminosity 
of 1.03 X lO'*" /i-~^ ergs "'^. These structures have features 
typical of radio relics IjGiovannini fc Ferretil [200J). They 
fit quite well on a projected ellipse with minor and major 
axis of ~ l.lh~^ and ~ 1.4/i~^ Mpc, respect ively. Adopt- 
ing a line-of-sight depth of ~ 189 h~^ kpc, as lBagchi et al.l 
(|200g), the three-dimensional ellipsoid has a volume V ~ 
0.15 h~^ Mpc^. Fro m the composi t e map of radio and X-ray 
emissions shown bv lBagchi et al.l (|2006l ) (see their fig. la), 
we infer that only ~ 20% of this volumen is filled by the 
radio relics (i.e. Kciic ~ 0.003 h~^ Mpc^). 

Some physical parameters of the cluster, relevant for 
the purpose of this work, are not provided by the obser- 
vations but can be obtained from simulations; in concrete, 
the gas density in the relics, nn, and the shock velocity, Vb- 
We consider a cosmological non-radiative hydrodynamical 
A''-body/SPH (Smoothed Particle Hydrodynamics) resimu- 
lation of a galaxy cluster that has been initially selected from 
a dark matter simulation for a standard ACD M cosmology 
with f^o = 0.3, h = 0.7, as = 0.9 and Qb = 0.04 (JDolag et all 
I2OO5I ). The presence of radio relics in Abell 3376 might be 
interpreted as the result of an on-going merger of subclus- 
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ters (|En/31in et al.lll998h . Thus, the chosen simulated clus- 
ter has similar dynamical state as the observed one; it has 

Mq. From the analy- 



a virial mass Mvir ~ 1.4 x 10 h 
sis of this simulation we obtain nn 



2 X 10" 



and 



Vb — 1000 km s~^. Considering that the temperature of the 

ICM, where the sh ock is propagating, is TIcm ~ 10~^rrciic 

Hoeft et al. 2004) the Mach number results 



M 



4.2 



Gabici fc Blasill2003l ). For the rest of the paper we express 
all numerical values adopting h = 0.7. 



3 CONTENT OF RELATIVISTIC PARTICLES 

Shock waves generated during the formation and evolu- 
tion of galaxy clusters are the main source for thermaliza- 
tion of the intraclust er gas and the acceleration of parti- 
cles ()Pfrommer et al.l 12006 1. The activity of radio galaxies 
embedded in the clusters also contribute to the population 
of relativistic particles, leaving fossil radio plasma that is 
detected as cavities i n X-ray surface brightness maps (e.g. 
IChurazov et al.l l200d ). Two radio galaxies have been ob- 
served in the cluster Abell 3376. The radio source MRC 
0600-399 is associated with the second brightest cluster 
member galaxy, and the other radio source is possibly orig- 
inated from an elliptical galaxy. These radio galaxies are 
located within the central region of the cluster from where 
thermal Bremsstrahlung emission is detected. Therefore, we 
can consider that the ring-shape radio structures present in 
the periphery of the cluster are not connected with these 
point sources. 

The morphology of the X-ray and radio emission ob- 
served in the cluster Abell 3376 suggests that it is un- 
dergoing a merg er. As numerical simulations show (e.g. 
iHoeft et al.ll2004l '). shock waves propagate in both directions 
along the line that connects the centres of the merging clus- 
ters, with the radio relics observed almost exclusively at the 
location of the shock fronts. Thus, particle acceleration in- 
duced by shock waves generated during this process is a 
suitable scenario for explaining the origin of the relics. 

Taking into account observational evidences and numer- 
ical results, we assume that the content of relativistic par- 
ticles in the relics arises as a result of the acceleration by 
mergers shocks, neglecting the possible contribution of the 
radio galaxies. In the following subsections, we describe the 
acceleration processes and losses that affect both electrons 
and protons, which determine the particle distributions and 
their subsequent evolution. 



3.1 Acceleration and losses 

The diffuse radio emission produced at the location of relics 
supports the presence of relativistic electrons and magnetic 
fields. Although there is no observational evidence of the ex- 
istence of relativistic protons in the radio relics, these par- 
ticles could be accelerated in the same way as electrons are. 
We focus here on the acceleration of particles at the 
(non-relativistic) shock front via a diffusive process, such as 
first order Fermi mechanism. In this theory, the time-scale 
related to the rate of energy gain of particles is 



7 



riE 
eBc 



where 7 is the Lorentz factor and E is the energy [E — 
'ymc?) up to which the particle, electron or proton, is 
accelerated; B is the magnetic field. The parameter r] = 
2-Kfsc{c/vsY involves the velocity of the shock, Vs, and 
the ratio /sc between the mean free path of the particle 
and its gyroradius; the latter can be expressed as r^ = 
S[eV]/(300 Z B[G\), being Z the charge number of the par- 
ticle. For the shock velocity we assume a typical value of 
Da ~ 1000 kms~^, as it is inferred from the analysis of our 
cluster simulation, being in agreement w ith general simu- 
lations results fe.g. IPfrommer et al.ll2006r ). Then, assuming 
Bohm diffusion (i.e. /ac = 1), we have ri = 5.7 x 10*. 

The rate of accelerated particles, Q(7), injected in the 
source follows a power-law energy distribution with a spec- 
tral index F = 2.1, typical of a diffusive acceleration mech- 
anism and in accordance with the estimated Mach number 
M (see Blasi 2003 for similar values of the spectral index in 
Coma cluster). Once the accelerated particles are injected, 
their spectral energy distributions evolve as a consequence of 
radiative and advection losses. The latter refer to the escape 
of particles from the acceleration region. 

In order to study the time-evolution of the energy dis- 
tribution of particles, we need to calculate the leptonic and 
hadronic radiative losses. In the following, we refer as pri- 
mary particles to those relativistic particles that have been 
accelerated at the shock front by the Fermi-like mechanism. 
Secondary particles are the products originated from the in- 
teractions suffered by relativistic protons with target ther- 
mal protons of the ICM. These inelastic proton-proton [pp) 
collisions are the main channel of hadronic losses. 

The cooling rate of relativistic protons as a result of pp 
interactions can be estimated by 



7pp = —4.5 X 10 



^nH [0.95 -I- 0.06 In ("-l- 



7 s 



(2) 



(jMannheim fc Schlickeiseij|l994h . where 7 is the Lorentz fac- 
tor for protons and nn is the average hydrogen density of 
the relics obtained from the simulation. Part of the energy 
lost by relativistic protons is used to create neutral pious 
which subsequently decay into 7-rays. The rest is used in the 
creation of charged pious that will decay finally producing 
electron-positron pairs (e ) and neutrinos. The secondary 
pairs are cooled by the same radiative processes that affect 
primary electrons. Lepton energy losses are estimated taking 
into account relativistic Bre msstrahlung, IC interactions a nd 
synchrotron radiation (e.g. ICinzburg fc Svrovatskiill 19641 ) . 

Since the ICM is fully ionized, we consider 
Bremsstrahlung losses for bared nuclei, which are given by 



7Br = -6.9 X 10 ^^nn^^ [ln7 -I- 0.36] 7 s"\ 



(3) 



Losses due to IC interactions depend on the energy den- 
sity of photons, f/ph. In the Thomson regime, they are given 
by 



71C 



-3.2 X 10~**[/ph7^ s"\ 



(4) 



(1) 



The two photon fields relevant for the physical processes 
taking place in clusters of galaxies are the CMB radiation 
and thermal X-ray photons originated within the ICM itself 
through thermal Bremsstrahlung. The former is given by 
t/cMB = aT"^ c~^ (1 -I- z)*, where a is the Stefan-Boltzmann 
constant and T = 2.725 K is the CMB temperature; thus, 
t/cMB — 1.2 X 10~^^ergcm~^ at the cluster redshift. The 
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latter is expressed as Ux ~ Lx (47r_R^ c)~^ (1 + z;)'', where 
Lx is the X-ray luminosity and 7? is the radius from which 
this luminosity is measured, which is about ^ 0. 15 — 0.3 
per cent of the virial radius (JBalestra et al.l 120071 ). In the 
case of the cluster under study, Lx — 2.44 x 10**ergs~^ 
and R « 0.3 Mpc, giving Ux « 1-3 x 10~^^ ergcm"^. Since 
Ux is two orders of magnitude lower than Ucmb, we can 
consider Ux negligible (in both Thomson and Klein-Nishina 
regimes of IC interactions) and Uph ~ Ucmb- Taking into 
account the energy for which the emission of the CMB spec- 
trum reaches its maximum (-Ecmb ~ 1.9 x 10"'* eV), leptons 
should reach energies greater than 1.6 x 10^^ eV in order to 
satisfy the condition ^ — EeEcMB / i'meC^)^ ^ 1, that defines 
the Klein-Nishina regime of the IC interactions. As we show 
in Section 1221 these high energies are never achieved by pri- 
mary electrons; hence, IC scattering developes in Thomson 
regime. 

Losses produced by synchrotron radiation are computed 



'Ysyncli 



-1.9 X W'^B^'Y^ s"\ 



(5) 



The calculation of synchrotron losses requires the knowledge 
of magnetic field, B. We estimate B by assuming equipar- 
tition between the energy density of the magnetic field and 
the relativistic particles, that is. 



B _ 

J-, — f^ei \ tip ~r '^62 5 
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where u^^ and Us2 are the energy densities of primary elec- 
trons and secondary pairs, respectively, and Up is the en- 
ergy density of protons. The energy density of i-particles 
(i = ei,p, 62) is expressed as 



E, n{Ei) dE,, 



(7) 



where n{Ei) is the particle energy distribution. For protons, 
n{Ep) = KpEp'^-^ exp{-Ep/E^'^'') but for leptons n^Ee^^^) 



is a broken power-law: 



, „ , 1 Ke E, 
"(^'=) = ^ KlE 
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(8) 



The spectral indices are Fej =2.1 for primary electrons 
and Fej = 2.08 for secondary pairs; E^, is the break en- 
ergy (see Section [3.31 for the estimation of Eh and K'^). In 
all cases _Ef""' is the maximum energy that particles can 
reach, that depends on the mechanism of energy gain and 
losses that affect the particles. The contribution of E^'^^ to 
integral ((TJ can be neglected because the slope of the par- 
ticle distributions is negative: n{Ei) ex E~ \ On the other 
hand, we assume a minimum energy -E™'" — 1 MeV for pri- 
mary electrons and iS™'" = 1 GeV for protons. Thus, for 
secondary pairs we estimated E™'" ~ Ep''^/6 ~ 0.17 GeV. 
Since the number of accelerated protons is unknown, we con- 
sider Up = aUei , where a is a free parameter in our model. 
We adopt three different values: a = (no proton acceler- 
ation), a — 1 (equal energy density in protons and in elec- 
trons), and a = 100 (protons dominate the energy density 
budget). 

Efficient production of secondary electrons requires a 
dense population of thermal protons as a target. However, 
this condition is not found at the peripheral cluster re- 
gions which are characterized by low density gas (nu ~ 



Table 1. Magnetic field, B, and normalization constant of the 
energy distribution of primary accelerated electrons, K^i , and 
protons, Kp, for the three different cases considered characterized 
by the parameter a. 
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9 X 10 
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[erg'"~^ cm""^! 



Kr, 



1 1.1 X lO""^ 

100 3.4 X 10-6 



1.4 X 10--^" 
1.1 X 10-15 
3 X IQ-i'^ 



[erg'" 



8.7 X 10-16 
2.5 X 10-14 
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Figure 1. Rate of energy gain of primary electrons for rj ~ 5.7 X 
10'^ compared with energy loss rates obtained for the case a = 1. 



10~^ cm"''). Therefore, the contribution of secondary pairs 
to the total energy density (Eq. [6]) is neglected. Then, ap- 
plying the constraint of the radio flux observed for Abell 
3376 (Fi, = 302 mjy), we estimate the magnetic field and 
the normalization constant of the energy distribution of pri- 
mary accelerated electrons, -ft'ei, and protons, Kp. Results 
are shown in Table [l] Using the estimated spectral energy 
distribution of protons, we obtain a ratio between energy 
densities iiea/up ~ 10~^, which justifies our previous as- 
sumption about the secondaries. 



3.2 Maximum energies 

The theoretical maximum energy of primary particles is de- 
termined by the competition between the rates of energy 
gain and losses, as described in the previous section. Fig. [T] 
shows the time-scales of energy losses of primary electrons as 
a function of their energy (Eqs.|31l3]and[SJ. These time-scales 
are compared with the one corresponding to the accelera- 
tion rate, given by Eq. [l] From this comparison, we obtain 
the maximum energy allowed for electrons, -E™^'', which is 
~ 9 X 10^^ eV. For the case of protons, we take into account 
the energy loss rate given by pp interactions (Eq. [Sjl , consid- 
ering that they have the same acceleration rate as electrons. 
We estimate proton maximum energies E™'"' ~ 1.8 x lO^'^ eV 
and 4.6 x 10^^ eV for cases with a — 1 and 100, respec- 
tively. At proton energies higher than ~ 5 x 10^* eV, pho- 
topair production against CMB photons starts to be im- 
portant, and beyond energie s of ~ 5 x 10^^ eV photopion 
losses become dominant (e .g. iBerezinskv fc Grigorevalll988l : 
iKelner fc Aharonianll2008l ). However, in the current analysis 
we do not need to take these losses into account since, as we 
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Table 2. Maximuin energies obtained for primary electrons and 
protons accelerated in the radio relics of Abell 3376 are given in 
the first and second columns, respectively. The maximum energy 
for secondary pairs e^ is presented in the last column. 



[cV] 



pmax 



9 X 10 



1 9.3 X IQi^ 
100 8.8 X 10^^ 



cV] 



/Ttmax 



cV] 



5.0 X lOi'' 4 X 10^^ 
1.3 X 10^** lOi'' 



will see below, the actual maximum energy of protons in the 
relics will be determined by additional constraints. 

The values of maximum energy found for both electrons 
and protons are valid as long as they allow the particles to 
remain within the acceleration region. Thus, the particles 
have to satisfy the constraint r^ < I, where r^ is the gy- 
roradius of the particle and I is the size of the acceleration 
region. By assuming that particles are accelerated by the 
shock waves traced by the observed radio relics, we adopt 
the approximation I ~ licUc ~ 0.3 Mpc; this value is ob- 
tained from the examination of the projected radio map of 
the cluster Abell 3376. For the electrons with energies close 
to the maximum value -E™^", we find that rg ~ 0.1 pc. Thus, 
the most energetic electrons are contained within the accel- 
eration region. The situation for protons is different since the 
corresponding maximum energy, E'™'"', makes their gyrora- 
dius larger than ZrcUc- Hence, iS™^" is determined from the 
condition rg = IrcUc- This estimation gives E^"^"^ ~ 9.7 x 10^'' 
eV and 2.6 x 10^° eV, for a = 1 and 100, respectively. How- 
ever, the time required for protons to reach these energies 
is higher than the lifetime of the relic (rj^ock ~ 1 Gyr; see 
the next subsection). Thus, the actual maximum energy to 
which protons can be accelerated is obtained equating race 
and Trciic. The resulting values of iJ™^" obtained from this 
constraint are lower, being ~ 5 x 10^^ eV and ~ 1.3 x 10^* 
eV for a = 1 and 100, respectively. 

The contribution of secondary pairs to the total spectral 
energy distribution (SED) will be taken into account even 
though the ambient conditions at the cluster outskirts are 
not favorable for the generation of these particles through 
pp interactions. To this aim, we first calculate the injected 
spectrum of secondary pairs taking into account the esti- 
mated distribution of relativistic protons, n{Ep), and the 
new parametrization of the inelastic cross-section of pp in- 
teractions 



o-inei(£'p) = 34.3 + 1.88L + 0.25L^ mb. 



(9) 



given by Kelner, Aharonian & Vugayov (2006), where L = 
ln(_E'p/l TeV). The maximum energies for secondary pairs 
are assumed to be equal to t he maximum energ y of the pho- 
tons produced by Tr^-decay (|Kelner et al.ll200d Fig. 12). 

We now analyse the time-evolution of the spectral en- 
ergy distribution of the three kind of relativistic particles 
present in the radio relics. 



3.3 Time-evolution of the energy distributions of 
relativistic particles 

Relativistic particles are injected in the radio emitting re- 
gion of the cluster with an energy distribution that follows a 
power-law, which subsequently evolves as a consequence of 
the radiative losses suffered by the particles. Such evolution 
is estimated, as a first order approximation, by adopting a 
model in which the energy distribution of the particles is 
determined by the following transport equation 

dn {t, 7) /dt + djn (t, 7) /d-y + n (t, 7) /tosc = Q(t, 7) (10) 



fe.g.. iKhangulvan et al■ll2007^ ■ where t is the time, Q{t,'y) 
is the particle injection function, wich for present case will 
be time independent, and Tcsc is the advection escape time 
defined as Tcsc ~ ^rcUc/'^adv In the shock rest-frame, the 
downstream (post-shock) plasma velocity is given by «adv ~ 
i'shock/4; thus, we obtain Tosc ~ 1.1 Gyr. 

Particles are injected into the relics with a rate Q{"/) as 
the shock wave propagates from the centre of the cluster up 
to the location of the relics. For the purpose of the present 
work, the variable t corresponds to different values of the 
lifetime of the shock, Tshock- When the shock reaches the lo- 
cation of the observed radio relics, this lifetime is referred to 
as Tg^ock- From the assumed shock velocity Vb ~ 1000 kms~^ 
(Subsection l3.1|l and the position of the observed radio relics 
in cluster Abell 3376, we have r!:^ack ~ 0.95 Gyr. This value 
is consistent with the one o btained from num erical simula- 
tions of cluster mergers (e.g. iHoeft et al.ll2004h . 

The time-derivative 7 is a function accounting for all the 
energy losses afi'ecting leptons and protons, that is, mainly 
synchrotron and IC interactions in the former case and pp 
interactions in the latter one. The solution of Eq. (|10|l for 
different shock lifetimes is 

7ctf 
7 

where 7off and t(7,7') are determined by 



'^shock 



d'y' 



and r(7, 7 ) 



d'y" 

7" I ' 



The cooling time of leptons is Tcooi = 7/710SS 



(12) 



2 x 



10 7~ s. If Tcooi > Tshock at a certain energy, the particle 
energy distribution has the shape characterized by Fe = F 
and Ke (see Table [l|. There is a break energy. Eh, above of 
which Tcooi < Tshock. lu this case, both primary electrons and 
secondary pairs suff'er radiative losses that affect the injected 
particle spectrum. As a consequence of this, the cooled spec- 
trum has an index F^^ ^ = Tei 2 + 1 with a normalization 
costant K'^^ ^ ~ EhKe^ 2 for energies Ee-^ 2 > ^b- The com- 
puted particle energy distributions of primary electrons and 
secondary pairs are shown in Fig. [2] As can be seen, the 
spectrum is broken at iJb ~ 5 x 10^ — 5 x 10^" eV for values 
of Tshock comprised in the range 100 Myr to 1 Gyr, being 
the break energies smaller for larger lifetimes of the shock. 
In this scenario, particles will reach the steady regime 
when Tshock > Tesc- For this reason, as rl^ada ~ Tesc, we can 
consider that by the time when the radio relic is observed, 
the particle energy distribution 71(7) is already steady. 
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Figure 2. Spectral energy distributions of primary electrons and secondary pairs for a = 1 (upper panels) and a = 100 (lower panels). 
Distributions for different particle injection time-scales (rghock) are shown in each panel. Different curves correspond to the following 
values of rsijtjck- 0.125 Gyr (dotted line), 0.25 Gyr (dot-dashed line), 0.5 Gyr (dashed line) and 1 Gyr (solid line). The steady regime is 
reached for ages ' 



relic 
shock * 



Regarding protons, the escape time-scale (tcsc ~ 
1.1 Gyr) and shock lifetime (rJJock ~ 1 Gyr) are shorter 
than the cooling time (rpp ~ 10'^ Gyr). Consequently, the 
energy distribution in the context of our scenario keeps the 
same spectral shape as that of the injected one. Since Tcsc is 
of the order of rJJock > "^6 can consider that the energy distri- 
bution of protons also reaches the steady regime at Tghock- 

All these calculations give us the steady energy distri- 
bution of the different population of relativistic particles, i.e. 
primary electrons, secondary pairs, and protons, that have 
become part of the relics. Electrons are responsible for the 
radio emission of the cluster Abell 3376 observed in its giant 
ringlike radio structures. As we have shown, this emission 
has been used to estimate the magnetic field strength in the 
acceleration region, and the normalization and maximum 
energies involved in the power-laws that represent the parti- 
cle energy distributions. In the following section, we estimate 
the contribution of both relativistic leptons and protons to 
the high energy emission of the cluster, evaluating the de- 
tectability of Abell 3376 in the MeV-TeV electromagnetic 
range. 



4 PRODUCTION OF GAMMA-RAYS AND 
LOWER ENERGY RADIATION 

Using the steady distributions of particles estimated in the 
previous section, we calculate the SEDs for the three dif- 



ferent combinations of proton and lepton relativistic energy 
densities (i.e. cases a = 0,1, 100). We take into account the 
most relevant non-thermal radiative processes according to 
the conditions of the ambient medium in the cluster Abell 
3376. 



4.1 Leptonic emission 

The differential emissivities q^{E.y) produced by leptons 
for synchrotron radiation, IC scattering and relativistic 
Bremsstrahlu ng are calculated usi ng the standard formu- 



lae gi ven by iBlumenthal fc Gouldl (119701 ) and IPacholczvkl 
(|l970l ). The luminosity produced in the relics is given by 
the equation 



E-yLSy = E^ q-.f{E-y) Vrolic, 



(13) 



where E^ is the photon energy and Vreiic is the volume of 
the emitting region. For the latter, we adopt a value of 
~ 0.09 Mpc^, as we have mentioned in Section (2] In the 
cases where a population of accelerated protons is assumed 
to be present (a = 1, 100), we estimate the contribution of 
secondary pairs to the SEDs for completitude. 

Results of the luminosity produced in the case where 
only primary electrons are accelerated, characterized by a = 
0, are shown in Fig.O We can see that IC interactions are the 
dominant process, with a luminosity Lie ~ 9.1x 10*^ ergs~^, 
at energies E^ Si 0.1 MeV and with a cut-off at ~ 10 TeV. 
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Figure 3. SED for a pure leptonic case (a = 0), represented 
by a solid line. The contribution of different radiative processes 
are shown: synchrotron radiation (dot-dashed line), IC scattering 
(dotted line) and relativistic Bremsstrahlung (dashed line). The 
VLA observation at 1.4 GHz is also included. 



Figure 4. SED for a mixed primary population of relativistic 
electrons and protons characterized by a = 1, represented by a 
solid line. The contribution of different radiation processes are 
identified as in Fig. |3] The contribution from secondary pairs is 
also included. 



The luminosity produced by relativistic Bremsstrahlung is 
negligible, in agreement with results shown in Fig. [T] 



4.2 Hadronic emission 

Proton-proton interactions produce 7-ray emission through 
TT^-decay. The corresponding differential 7-ray emissivity is 
calculated as 



liiEj) 



qno{En) 



"(Bt) 



vw 



dE^, 



(14) 



where .B™'"(.B^) ^ E^ + mlc^/AE-,. Applying the 5- 
functional approximation fo r the differential cross sectiorij 
IjAharonian fc Atovanll200Cl ). the pion emissivity becomes 



q^o{E„) ^^nuJp 



[m^c^ + S^) 



{m,c^ + ^) (15) 



for proton energies greater than the energy threshold Eth ~ 
1.22 GeV and lower than 100 GeV. Here, k is the mean 
fraction of the kinetic energy Skin = Ep — nipC^ of the proton 
transferred to a secondary meson per collision. For a broad 
energy region (GeV to TeV) we have that k ~ 0.17. The 
proton flux is given by Jp{Ep) = (47r/c) n(i5p). The total 
cross section of the inelastic pp collisions is giv en by Eg. Q. 
For E p > 100 GeV, the equations given by iKelner et al.l 
l|2006l ) are used. The specific luminosity is estimated in the 
same way as for electrons fEg. llSp . 

The results of our calculations for the case a — 1 are 
shown in Fig. |4] Just as the case charaterized by a = 0, 
the SED is dominated by the IC interactions, with a lumi- 
nosity Lie ~ 7.4 X 10*^ erg s~^. The emission at energies 
higher than ~ 1 GeV is produced by neutral pion decay, 
reaching a luminosity Lpp ~ 1.6 x 10"^* erg s~^ with a cut- 



off at E^ 



-|17 



10 eV. However, its contribution to the SED 



becomes evident at energies > 10 TeV. 



^ This approximation considers only the most energetic neutral 
pion that is produced in the pp reaction aside of a fireball com- 
posed by a certain num ber of less energetic 7 r-mes ons of each 
flavor. See a discussion in lPfrommer &: En/31inl ll2004f) . 



In the case corresponding to a = 100, the relativistic 
energy density of protons is higher than that of the pri- 
mary electrons. However, the proton dominance does not 
imply necessarily a significant increase of photon produc- 
tion from pp-interactions because of the low density of ther- 
mal protons at the location of radio relics. Hence, the power 
produced by Tr^-decay does not dominate the SEDs, as can 
be seen in Fig. [5] The luminosity Lpp is of the order of 
~ 4.2 X 10^^ ergs~^, being slightly larger than the one cor- 
responding to the case with a = 1. In contrast to the cases 
with a = and a = 1, the IC emission is lower than the 
synchrotron one, with luminosities Lie ~ 7.1 x 10^" ergs~^ 
and Lsynch ~ 3.8 x 10''^ergs~^, respectively. The fact that 
radiation produced by IC scattering is reduced for the case 
a = 100 might be explained as a consequence of equipar- 
tition between the magnetic field and the particle energy 
density. As can be seen in Table [l] the magnetic field be- 
comes larger as the value of a increases. This will reduce the 
amount of electron energy released via IC scattering with 
respect to that released via synchrotron. In addition, the 
available electron energy is also reduced to explain the radio 
fiuxes due to the larger magnetic field. All this makes the 
IC luminosity smaller. Concerning the emission produced 
by secondary pairs, we find that their contributions to the 
SEDs by different radiation processes are larger than that 
for the case a = 1, but are still irrelevant. This is consistent 
with the low efficiency of pp-interactions taking place in the 
outskirts of the cluster. 



4.3 Gamma-ray emission 

The SED shown in Figs. [3] to [S] are built considering the 
contribution of different non-thermal radiative processes in 
a wide range of energies, from radio to gamma-rays. In 
the present work, we are particularly interested in the de- 
tectability of the cluster Abell 3376 at gamma-ray energies. 
We find that the main processes contributing to the emission 
at these high energies are IC scattering and 7r°-decay, being 
the luminosity of the latter considerably smaller even for the 
case where a large energy density of protons is considered. 
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Figure 5. Idem as Fig|4] but for a case dominated by relativistic 
protons [a = 100). 



The HESS array could detect the estimated emission from 
A3376 in exposure times between one day to one month, for 
different values of a. On the other hand, more than one year 
is necessary to detect with GLAST the non-thermal gamma- 
ray emission estimated for the cluster studied in this work. 
However, small changes in the assumptions of our model can 
yield higher fluxes that would result in a considerably reduc- 
tion of the exposure time for detection. For instance, if sub- 
partition between magnetic and relativistic particle energy 
densities is considered instead of equipartition, there would 
be an increase of the IC emission. Thus, if the IC luminos- 
ity of Abell 3376 is Lie ~ 10*^ erg s~^ at energies above 
100 MeV it would be detected by GLAST in less than one 
year. (i.e. within the initial all-sky survey). Finally, photons 
produced by pp-interactions could be might by the future 
Cherenkov telescope HESS II on short timescales. 



5 SUMMARY AND PROSPECTS 

We have presented a model for the emission produced by 
the nearby cluster Abell 3376 in a wide range of ener- 
gies, from radio to gamma-rays, considering the contribu- 
tion of different non-thermal radiative processes. This clus- 
ter presents strong evidence of merger activity, characterized 
by large-scale ring-shaped synchr otron radio structu res, that 
are identified with radio relics (IBagchi et al.ll2006l ). These 
kind of radio structures detected in several clusters of galax- 
ies suggest a rich content of accelerated particles in the ICM, 
with energies as high as ~ 10^^ and ~ 10^* eV for elec- 
trons and protons, respectively. Taking into account that 
radio relics are tracers of merger shocks, we assume that a 
diffusive shock acceleration mechanism acting at the wave 
shocks gives place to such a population of relativistic par- 
ticles. These particles are subsequently cooled by different 
leptonic and hadronic radiative processes. 

We have considered syncrotron radiation, IC scatter- 
ing and relativistic Bremsstrahlung that affect relativistic 
primary electrons. Relativistic protons are involved in pp- 
interactions with the thermal protons of the ICM. This pro- 
cess give raise to a population of secondary pairs, e , that 
are cooled by the same processes that affect the population 
of primary electrons. The parameters involved in these non- 



thermal processes, such as the magnetic field and normal- 
ization of the energy distribution of accelerated particles, 
have been estimated from the constraints established by the 
observed radio power in the relics and the assumption of 
equipartition. X-ray observations allows us to extract infor- 
mation about the volume of these acceleration sources. On 
the other hand, results from a simulated cluster have pro- 
vided the typical gas density of the relics and the velocity 
of ths shock. 

We find that particle energy distributions reach a sta- 
tionary regime by the lifetime of the relics (~ 1 Gyr). The 
steady energy distribution of particles are computed and 
the SEDs of the radiation produced at the relic position de- 
rived. In our calculations, we have adopted different ratios 
(a = 0, 1, 100) for the relativistic proton/electron energy 
densities. At larger energies, gamma-ray emission is mainly 
produced by vr^-decay but with a considerably smaller lu- 
minosity, even for the case where a large energy density of 
protons is considered [Lpp ~ 4 x 10^^ erg s~^). This is con- 
sistent with the fact that no cluster has been yet observed 
at these high energies. 

The proximity of the cluster studied in this paper and 
its high content of relativistic particles in the radio-relics 
make Abell 3376 an interesting potential target for the in- 
vestigation of gamma-ray emission in this type of objects. 
This source might be detected at 7-rays by GLAST satel- 
lite and by HESS array with a reasonable exposure time. In 
addition, future Cherenkov telescopes in the southern hemi- 
sphere, like HESS II, could easily detect the high-energy 
emission from Abell 3376. 
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